Critical dynamics of singlet and triplet excitations in strongly frustrated spin systems 
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New data for the two-dimensional J\ — Ji - model shows that the critical points for singlet and 
triplet excitations near J2/J1 ~ 0.38 are very close together, or coincident. We propose that this 
is a more general result: in strongly frustrated spin systems the critical dynamics of singlet and 
triplet excitations are closely related. An effective field theory for such a critical point, or points, is 
developed, and its implications discussed. The analysis suggests the existence of gapless or almost 
gapless singlet excitations in the magnetically disordered phase. 



PACS: 75.30.Kz, 75.30.Ds, 75.fO.Jm, 75.40.Gb 

Quantum phase transitions in quantum antiferromag- 
nets have attracted a great deal of attention during the 
last decade. This field contains numerous previously un- 
known phenomena. Investigation of quantum phase tran- 
sitions in quantum antiferromagnets is very important 
for quantum magnetism, but is of even wider importance 
because it allows one to construct and to analyze effec- 
tive nonperturbative field theories that can be used in 
other parts of physics. In this paper we concentrate on 
rotationally invariant (2+l)D systems. Transitions de- 
scribed by the nonlinear O(N) field theory are very well 
understood, see, e.g., Refs. |l]-|3| The most widely known 
example is the 0(3) theory that describes a real vector 
field v with Lagrangian density 
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The effective mass vanishes at some critical point, to 2 = 
a v{g — gcv), where g is some generalized "coupling con- 
stant", and a v is some coefficient. It has been argued 
by Haldane 0] and Chakravarty, Halperin and Nelson 
p| that this theory describes a quantum phase transition 
between the Neel state (g < g cv ) and the magnetically 
disordered state (g > g cv ) in a two-dimensional quantum 
antiferromagnet. Perhaps the simplest concrete example 
of such a system is the Heisenberg antiferromagnet on 
a double layer square lattice. In this case the "coupling 
constant" is equal to the ratio of interlayer and inlayer 
antiferromagnetic interactions, g = Jj_/J, g cv ~ 2.525, 
see Ref. Jp|. For g > g cv the excitation spectrum of the 
system (|l|) is gapped, and for g < g cv the spectrum is 
gapless because of Goldstone spin waves. 

The 0(1) field theory describes a real scalar field s with 
Lagrangian 
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where m 2 s = a s (g cs — g). The field expectation value 
< s > is zero at g < g cs and nonzero at g > g cs . The 
spectrum is gapped in both cases. A simple example 
of such a system is a Heisenberg antiferromagnet on a 
"planar pyrochlore" lattice where the entire low energy 
dynamics is driven by spin singlet excitations M. 



The situation in strongly frustrated spin systems is 
much more complicated because the dynamics of triplet 
and singlet excitations are not independent. The 2D 
Ji — J2 model is an ideal testing ground for understanding 
of such dynamics. The Hamiltonian of the model reads: 



H = J\ Si • Sj + J2 /J Sj • Sj, 



(3) 



where J\ is the nearest-neighbor, and J2 is the frus- 
trating next-nearest-neighbor Heisenberg exchange on a 
square lattice. Both couplings are antiferromagnetic, i.e. 
Ji,2 > and the spins Si = 1/2. At small J2 the system 
exhibits Neel ordering and at small J\ it exhibits collinear 
ordering. There is also a general consensus on the disap- 
pearance of the magnetic ordering at 0.38 < J2/J1 < 0.6, 
see, e.g., Ref. ||. The structure of the magnetically dis- 
ordered state is of great interest. An important insight 
into the disordered regime was achieved by field-theory 
methods [p|,p| and dimer series expansions pOpdll]. The 



above works have shown that the ground state in this 
regime is dominated by short-range singlet (dimer) for- 
mation in a columnar pattern. More recent studies based 
on the Green function Monte Carlo method favored ei- 
ther a plaquette state or a dimer state with plaquette 
modulation [fl3"[ . And finally the same group that claimed 
stability of the plaquette state now claims stability 
of the translationally invariant RVB state Q|. The above 
series expansions JTo|,pd| j8| , as well as Monte Carlo numer- 
ical calculations ]12|-|14| in various forms, have difficulties 
in distinguishing between states with very close energies, 
as occurs in the J\ — J2 model. In such approaches one 
starts from some trial wave function (dimer, plaquette, 
RVB,...) and such an ansatz is always very restrictive: it 
is impossible to get a translationally invariant state us- 
ing the dimer or plaquette expansions, and on the other 
hand it is impossible to get a spontaneous violation of 
the translational symmetry starting from the trial RVB 
state that is explicitly translationally invariant. 

A novel approach was implicitly formulated in our work 
p"5[ , see also Ref. Q . The idea of the approach is to use 
numerical calculations only in the well established Neel 
phase at J2/J1 < 0.38. This calculation allow us to find 
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all the relevant quasiparticles in the system, and hence to 
determine the effective field theory describing the dynam- 
ics of the system. Since the transition is of the second 
order the same field theory describes the magnetically 
disordered phase and we have no problem to describe this 
phase. In the present work we explicitly formulate the 
approach, present new, much more accurate, numerical 
data, derive the effective Lagrangian of the model, and 
finally discuss implications of this effective Lagrangian. 

For the J\ — J2 model we define the generalized cou- 
pling constant as g = Jij J\. It is well established that 
at g < g cv there is Neel ordering in the system. There 
are 3 independent estimates for g cv . The first comes 
from dimer series expansions, g cv w 0.38, see Ref. [||. 
However this estimate is problematic since it already as- 
sumes spin dimer order at g > g cv . The second estimate, 
g cv = 0.4 ± 0.02, comes from Ising series expansions for 
the staggered magnetization in the Neel phase Jll| . The 
third estimate g cv — 0.38 ± 0.03 comes from the Ising se- 
ries expansions for spin-spin correlators in the Neel phase 
Using the second and the third estimates, and av- 
eraging them as independent ones we find 
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0.39 ±0.016^ 0.39 ±0.02. 
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The Neel ordering corresponds to the condensation of a 
vector field, therefore there is no doubt that some ef- 
fective Lagrangian determines the dynamics of the field 
both at g < g cv and g > g cv . 

However the vector excitation is not the only relevant 
degree of freedom in the system. There are also the 
low energy singlet excitations. We do not have a re- 
liable numerical technique for direct calculation of sin- 
glet excitations, but we do have a well developed se- 
ries expansion technique for calculation of static sus- 
ceptibilities. A static susceptibility is proportional to 
the corresponding Green's function at zero frequency 
Xq Gq(uj = 0) ~ Z q /Wq, where cjq is the quasiparti- 
cle energy, and Z q is the quasiparticle residue. So at the 
critical point 1/x must vanish approximately as (g — g c ) 7 , 
with 7 = v(2 — rj), where v is the critical index for the 
spectral gap, A oc {g — g c ) v , and vr\ is the critical index 
for the quasiparticle residue, Z oc (g — g c ) UTI ■ For the 
simple 0(1) Lagrangian (|J) the critical index is 7 w 1.2 
(see, e.g., Ref. g]). 

In this paper we use Ising series expansion fll|| tech- 
niques to compute the susceptibility of this system. The 
linked-cluster expansion method has been reviewed in a 
recent article |fL7 |, and will not be repeated here. 

In Ref. H we have calculated the susceptibility of the 
Neel state with respect to the external field. 
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where the coefficient / gives the strength of the field. The 
susceptibility is defined in the usual way 
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where E is the ground state energy. This susceptibility 
probes spin singlet excitations that correspond to a spin 
columnar dimerization. The susceptibility is diverging 
at g = J2/J1 ~ 0.4 and this is a clear indication of the 
critical singlet excitation. 

Unfortunately convergence of the series expansion for 
X is somewhat erratic and therefore the error bars in the 
numerical data [|l5| are rather large. In the present work 
we calculate susceptibilities Xz and x± with respect to 
the fields 

F dz = fJ2(-^S%S^ 3l (7) 

Here z is direction of the staggered magnetization in the 
Neel state. Both Xz and x± probe the same excitations 
as x- Using the Ising expansion about the Neel order 
JlTf , we have computed the series for Xz and x± in pow- 
ers of exchange anisotropy to order 7 for given values of 
J2/J1. The calculations involve a list of 28474 distinct 
connected clusters up to 9 sites. We do not display the 
series here but can provide them on request. Fortunately 
convergence of the series for x± and Xz is much better 
than that for %. Let us note that the field Fd can be rep- 
resented as Fd = Fd z + Fd± . However the susceptibility 
X cannot be represented as a simple superposition of Xz 
and x±i X Xz+X±- Therefore a relatively poor conver- 
gence of the series for % does not contradict the better 
convergence for Xz and The results for 1/xz and 

1/X-L. are presented in Fig.l and Fig. 2. The susceptibili- 
ties are diverging at the critical point g cs . From Fig.l we 
conclude that g cs = 0.38 ± 0.04, and from Fig. 2 we con- 
clude that g cs = 0.41 ±0.03. Considering these estimates 
as independent ones and averaging them we find 
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= 0.40 ± 0.024 -> 0.40 ± 0.02. 
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We have also calculated the susceptibility Xp with respect 
to the field 

F P = EM)^ (S%S$ W + S%S% W ) , (9) 

that probes excitations of plaquette symmetry, but this 
is not coupled to the spin-dimer excitations. The Ising 
series for Xp has been computed up to order 7. The calcu- 
lations involve a list of connected clusters different from 
that for Xz and x± because the symmetry of the field 
is different. The number of cluster in this calculation is 
28450, also up to 9 sites. The results of the Ising series 
expansion calculations for 1/Xp are shown in Fig. 3. Ac- 
cording to the data the plaquette excitation is not critical 
and hence is not relevant. 
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The susceptibilities (Q) probe the real scalar field s 
that corresponds to the spin columnar dimerization. The 
ground state expectation value of the field is zero at 
g < g cs (no spontaneous dimerization), but the suscep- 
tibilities are diverging at g cs . Therefore one can say the 
effective Lagrangian for the field s is given by Eq. <^). 
This is a possible scenario, but does not explain why the 
critical points for the scalar field (||) and for the vector 
field (|4|) coincide or at least are very close: in this scenario 
this is just a pure coincidence. To be absolutely correct 
we have to say that there is another scalar field s' that 
corresponds to the replacement i — ► j, j — > i in Eqs. ([?]). 
So the symmetry of the transition is 0(1) x O(l). How- 
ever this additional s'-fiels is completely identical to s, 
and therefore it does not require a special consideration. 

The fields v, s and s' are not independent: they cer- 
tainly interact with each other. Therefore there is an al- 
ternative scenario that naturally explains the coincidence 
of the quantum critical points. Consider the effective La- 
grangian of an interacting massive vector field and two 
massless scalar fields. 

£ = \ [(<9 t v) 2 - c 2 (Vv) 2 - m 2 v 2 ] - ^ (v 2 ) 2 (10) 
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where 



a v (g — g C v)- This Lagrangian incorporates 



all the relevant degrees of freedom. At g < g cv the vector 
field develops the ground state expectation value < v >= 
V 'dv(gcv — g)/u v . and, hence, because of the interaction 
the massless scalar fields attain an effective mass m 2 ff = 
d v (gcv — g)u sv I (2u v ) . Thus in this scenario the scalar 
susceptibilities are diverging exactly at the same point 
where the Neel order disappears, \ oc l/(g C v — <?)■ In 
the magnetically disordered phase at g > g cv both the 
scalar fields remain massless. This scenario corresponds 
to a translationally invariant spin liquid state. One can 
also introduce in Eq. ( |l0|) a small negative bare mass for 
the scalar fields, m 2 = — |m 2 |. However due to Eqs. (Q) 
and (^|) this mass must be very small because it splits the 
critical points. In this case there is a region between the 
critical points where the spin-dimer and the Neel order 
coexist H|J|. 

We mention, in this connection, that in one- 
dimensional systems the Lieb-Schultz-Mattis (LSM) the- 
orem guarantees that a uniform phase in a quantum spin 
system is gapless and a gapped phase breaks translational 
symmetry [Q. Extension of the LSM theorem to two di- 
mensions was proposed [ fl9[ but not proven (see also Ref. 
p0|). We note that the effective Lagrangian ( [Tof ) is in 
perfect agreement with the conjectured two-dimensional 
extension of the LSM theorem. If the bare singlet mass 
m 2 is zero then the magnetically disordered state is uni- 
form and there are two gapless excitations in the state. If 



to 2 = — |m 2 | ^ then one of the scalar fields develops the 
ground state expectation value < s >= \f\m1\/u s (spon- 
taneous dimerization) and another attains finite positive 

mass due to the interaction terms, m 2 f f = ^-\m 2 A. 

The effective Lagrangian (|10| ) corresponds to the trans- 
lationally invariant intermediate spin liquid state (SL) 
and in this part it is in agreement with the conclusion of 
the recent work However, as we have shown, this 

SL state must have gapless singlet excitations. There- 
fore: 1) the dimer-dimer correlation function must fall 
off as oc 1/r, where r is the separation between dimers, 
2) for a finite cluster of N sites the static susceptibility 
with respect to the dimerization field (^) must scale as 
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Both these conclusions disagree with results of Rcf . [Q . 
In our opinion this can be related to the variational na- 
ture of that calculation. We would like to stress that for 
calculation of the susceptibility the external field must 
be weak, fv~N E s — E ~ l/s/N. If this criterion is 
violated, the ground state energy shift is proportional to 
oc fV~N and the second derivative (|^) vanishes. 

It is more likely that the bare singlet mass m 2 — ~ |m 2 | 
is small, but non-zero. The dimer series expansions Q 
clearly show two nonzero dimer order parameters in the 
SL phase. To avoid misunderstanding we stress that we 
do not claim that one can derive a translationally invari- 
ant ground state starting from the dimer basis. This is 
certainly impossible. However if the true state is a trans- 
lationally invariant one then the dimer series must gen- 
erate zero dimer order parameters as an indication that 
there is something wrong with the dimer state. This does 
not happen. In principle it is certainly possible that con- 
vergence of series is so peculiar that the standard Pade 
approximant method based on a limited number of terms 
does not give the right answer, but this is an unlikely 
event and there no examples of such behavior. In the 
case of the dimer ground state the static susceptibility of 
a finite cluster with respect to the dimerization field (|5|) 
must scale as 
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For the calculation the external field must be very weak, 
fN <C AE ~ exp(— aN). If this criterion is violated, the 
ground state energy shift is proportional to oc fN and 
the second derivative (^J) vanishes. 

In conclusion we have demonstrated that the critical 
dynamics of singlet and triplet excitations in the two- 
dimensional J\ — J2-model are closely related. The criti- 
cal points for singlet and triplet excitations are very close 
or coincide. This is a new type of critical behavior and 



3 



we argue that this is a generic property of strongly frus- 
trated spin systems. We have formulated the effective 
field theory describing the critical behavior. The analysis 
favors gapless or almost gapless singlet excitations in the 
magnetically disordered phase. The effective field theory 
deduced in our work supports the two-dimensional ex- 
tension of the Lieb-Schultz-Mattis theorem conjectured 
some time ago. 
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Advanced Computing (APAC) National Facility. 




FIG. 1. The plot of 1/Xzi where \z ts the dimer suscepti- 
bility calculated in the Neel state using Ising series expansion. 
The value of 1/xz vanishes at g cs = 0.38 ± 0.04. 
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FIG. 2. The plot of 1/X-L, where x± is the dimer suscepti- 
bility calculated in the Neel state using Ising series expansion. 
The value of l/x± vanishes at g cs = 0.41 ± 0.03. 




FIG. 3. The plot ofl/xp, where \ p is the plaquette suscep- 
tibility calculated in the Neel state using Ising series expan- 
sion. The value of 1/xp does not vanish, so excitations with 
plaquette symmetry are not critical. 
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